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Multiple  stochastic  integral  expansions  are  applied  to  the  problem  of 
filtering  a  signal  observed  in  additive  noise*  It  is  shown  that  the  optimal 
mean- square  estimate  may  be  represented  as  a  ratio  of  two  multiple  integral 
series*  A  formula  for  expanding  the  product  of  two  multiple  integrals  is 
developed  and  applied  to  deriving  equations  for  the  kernels  of  best,  finite 
expansion  approximations  to  the  optimal  filter.  These  equations  are  studied 
in  detail  in  the  quadratic  case. 
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SIGNIFICANCE  AND  EXPLANATION 


A  common  problem  in  the  analysis  of  stochastic  systems  is  the  estimation 
of  a  stochastic  process  given  only  noise-corrupted  or  incomplete  observations. 
Examples  occur  in  communications  theory  when  one  wants  to  estimate  a  signal 
sent  over  a  noisy  channel  or  in  time  series  problems.  If  x(t)  is  a 
stochastic  process  denoting  the  signal,  the  observations  are  typically 
modelled  by 

y(t)  =  jJ  h(x(s))ds  +  dW(t)  , 

where  W(t)  is  an  independent  increments  "noise"  process,  usually  Brownian 
motion.  The  problem  of  filtering  is  to  build  an  estimate,  i.e.,  filter,  of 
x(t)  using  the  observations  y(s),  s  <  t.  Theoretical  characterizations  of 
best  mean-square  estimates  are  known,  but  can  be  translated  into  effective 
solutions  only  in  special  instances.  In  this  paper,  the  general  filtering 
problem  is  treated  by  attempting  to  expand  filters  in  series  of  multiple 
stochastic  integrals  of  the  form 

Jj  /o1— C  a(t'S1 . sr)dy(sr)...dy(s1)  . 

Two  primary  issues  raised  by  this  idea  are  considered;  representation  of  the 
optimal  mean-square  estimate  by  multiple  integral  expansions,  and  construction 
of  suboptimal  estimates  using  a  finite  number  of  multiple  integrals.  It  is 
shown  that  expansion  of  the  optimal  filter  is  indeed  possible,  and  a  method  is 
presented  for  finding  best,  finite  expansion  estimates.  A  rudimentary  algebra 
of  multiple  integral  expansions  is  first  developed  as  a  tool  to  prove  these 
results. 

The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


MULTIPLE  INTEGRAL  EXPANSIONS  FOR  NONLINEAR  FILTERING* 

Daniel  Ocone 

t. 1  Introduction 

In  the  additive  noise  model  of  filtering,  information  about  a  stochastic  process 
x(t),  t  >  0,  called  the  signal,  is  received  through  observations  of  the  form 

y(t)  =*  h(x(s))ds  +  w(t)  t  >  0 

w(t)  is  a  noise  term  that  corrupts  the  signal,  and  it  is  usually  assumed  to  be  a  Brownian 
motion.  The  filtering  problem  is  to  estimate  from  the  observations  y(s),  0  <  s  <  t,  a 
given  moment  f(x(t>)  of  the  signal  at  time  t,  and,  if  estimators  minimising  mean- 
square-error  are  desired,  this  means  calculating  the  conditional  mean  F{f(x(t))  IF  *}, 

F^  :  -  a{y(s)  |  0  <  s  <  t).  E{f(x(t))  |  F^}  is  henceforth  referred  to  as  the  optimal 

filter.  Two  fundamental  characterizations  of  the  optimal  filter  are  available:  a)  a  Bayes 
formula  for  E(f(x(t))  I  F^}  as  the  ratio  of  two  conditional,  functional  integrals 
(Kallianpur,  Striebel  [9),  cf.  $1.2  of  this  paper):  b) ,  in  the  case  that  x(t)  is  Markov, 
a  representation  of  the  optimal  filter  as  a  stochastic  integral  against  the  Innovations 
process,  v(t)  “  y(  t)  -  /j j  E{h(x(s))  I  F^}ds,  the  stochastic  integrand  being  adapted  to 
the  observation  process  (Fujisaki,  Kallianpur,  and  Kunita  [2]).  However,  though 
theoretically  deep,  these  results  lead  to  explicit  and  analytically  computable  solutions 
only  in  special  Instances. 


• 

This  work  formed  part  of  the  author'a  Ph.D.  thesis  in  Applied  Mathematics  at  M.I.T.  under 
the  supervision  of  Professor  Sanjoy  K.  Mltter. 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041.  This  material  is 
based  upon  work  supported  by  the  National  Science  Foundation  under  Grant  No.  MCS-7927062, 
by  a  National  Science  Foundation  Graduate  Fellowship,  and  by  the  Air  Force  Office  of 
Scientific  Research  under  Grant  No.  AFOSR-77-3231B. 


This  paper  studies  the  application  of  multiple  stochastic  integral  expansions  to  the 


filtering  problem.  Any  filter,  optimal  or  suboptimal  is  actually  an  anticipating 
functional  of  the  observation  process ,  thus  suggesting  that  filters  be  represented  and 
analyzed  within  a  framework  for  functional  expansions.  Multiple  stochastic  integrals  prove 
useful  for  this  purpose.  In  fact,  their  definition  originates  in  Wiener's  homogeneous 
chaos  theory,  which  constructs  orthogonal  decompositions  of  spaces  of  f inite-variance 
functionals  of  Gaussian  processes  (cf.  Kallianpur  [8]  and  Kida  [5]).  In  the  Brownian 
motion  case,  each  subspace  of  the  decomposition  corresponds  to  the  space  of  multiple 
stochastic  integrals  of  a  given  order,  and,  thus,  Wiener's  theory  shows  that  any  L  - 
functional  of  the  Brownian  motion  may  be  expanded  in  a  series  of  multiple  integrals. 
Multiple  integrals  have  been  used  already  to  solve  a  number  of  specific  estimation 
problems.  Marcus,  Mitter,  and  Ocone  (13]  apply  the  homogeneous  chaos  theory  to  compute 
conditional  statistics  of  polynomial  functionals  of  a  Gauss-Markov  process  observed  in 
white  noise,  and  Hida  and  Kallianour  [6]  use  multiple  integrals  to  predict  non-linear 
functions  of  Brownian  signals  given  perfect  observations.  In  cumulant  approximations  of 
the  conditional  density  in  filtering,  Gterno  (1]  also  derives  expressions  using  multiple 
integrals.  Here,  we  seek  to  apply  multiple  integrals  of  the  form 


att.Sj 


,s  )dy(s  )»**dy(s  )  , 

n  n  i 


where  a(...)  is  deterministic,  to  the  general  filtering  problem.  We  focus  on  two  basic 
issues;  the  expansion  of  the  optimal  filter  by  expressions  involving  multiple  integrals, 
and  the  construction  of  best  suboptimal  filters  having  a  finite  multiple  integral  expansion 
of  specified  order.  It  is  important  to  observe  that  the  stochastic  integrals  we  employ  are 
formed  from  the  observation  process  and  not  the  innovations  process.  At  first.  Integration 
against  innovations  might  appear  to  be  an  attractive  idea  because  the  innovations  process 
is  Brownian,  integrals  of  different  orders  are  thus  orthogonal,  and  homogeneous  chaos 
theory  can  be  applied.  However,  in  practice  the  innovations  process  is  not  available  since 
its  construction  requires  the  estimate  E{h(x(t))  I  F^},  to  compute  which  is  generally  a 
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difficult  filtering  problem  itself.  Integrals  using  y(»)  directly  are  thus  more  natural, 
but,  due  to  their  more  general,  usually  non-Gaussian  character  are  more  difficult  to 
apply.  For  example,  in  suboptimal  estimation  one  might  like  to  project  random  variables  on 
a  sum  of  spaces  of  multiple  integrals.  This  is  easily  done  for  Brownian  Integrals,  using 
the  orthogonality  of  different  order  integrals  and  explicit  formulae  to  calculate  the 
integrands,  but  not  so  easily  for  more  general  integrals,  where  the  orthogonality  structure 
and  kernel  formulae  are  lost.  In  this  paper  we  describe  a  method  for  analyzing  y( •)  - 
based  integrals,  that,  in  particular,  allows  resolution  of  this  projection  problem. 

The  paper  is  organized  as  follows.  (1.2  Introduces  the  precise  filtering  model  we 
consider  and  recalls  the  Kallianpur-Striebel  formula  for  the  optimal  estimate.  X  central 
feature  of  this  formula  is  the  fact  that  the  y ( •)  process  is  absolutely  continuous  with 
respect  to  Brownian  motion.  Transformations  of  measure  so  that  y( •)  becomes  Brownian 
will  be  an  underlying  component  of  our  analysis  of  y( • ) -based  Integrals.  (2  is  a  self- 
contained  treatment  of  multiple  integrals  of  Brownian  and  observation  processes.  We  define 
multiple  stochastic  Integrals,  prove  technical  lemmas  for  later  use,  and  develop  some 
useful  properties  of  the  integrals.  Of  particular  importance  is  the  multiplication  formula 
<  theorem  2.1),  which  shows  how  to  express  the  product  of  multiple  Integrals  in  a  multiple 
integral  expansion,  thus  providing  a  rudimentary  algebra  for  handling  expansions.  We 
present  the  applications  to  filtering  in  section  3.  In  (3.1,  we  show  that  the  optimal 
filter  can  be  represented  as  the  ratio  of  two  multiple  integral  expansions,  essentially  by 
expanding  the  Kallianpur-Striebel  formula.  (3.2  addresses  the  issue  of  finding  the  best 
(mean  square)  estimate  of  the  form 


ao(t)  +  fo  +,”+  /o  /q1  ir  Wr***r>dy<V**‘dy(V 


By  combining  the  expansions  of  (3.1  and  the  multiplication  formula,  we  derive  a  system  of 
linear  integral  equations  for  the  kernels  {a^}^_0.  In  the  method  of  analysis  is 
to  transform  measures  to  a  space  on  which  y( •)  is  a  Brownian  process  and  then  to  apply 


the  multiplication  formula  to  discover  the  effect  of  the  Radon-Nikodym  derivative  so 
introduced.  The  remaining  sections  apply  these  results,  first  to  rederiving  the  Kalman 
filter,  second  to  finding  best  quadratic  filters. 

It  is  a  pleasure  to  thank  Professor  S.  K.  Hitter,  for  suggesting  this  problem  and  for 
inspiring  and  guiding  the  research. 

1.2  Filtering  preliminaries 

The  precise  filtering  model  to  be  considered  is  as  follows,  bet  the  underlying 
probability  space  be  denoted  P).  For  0  <  T  <  »,  let  {x(t)  I  t  e  fO,T)  }  be 

a  measurable,  real-valued  process  on  (n,F,P),  h(s,x)  a  Borel  function  on  [0,T]  x  R, 
and  w(t)  a  standard  Brownian  motion  on  ( ft,  ,P),  such  that 


Set 


i)  w<  •)  is  independent  of  x(  •) 

ii)  K  J*  h2(s,x(b))ds  <  «•  . 

'  o 


y(t)  -  /*  h(  S,x(  s) )  ds  +  w(t)  t  6  [0,T]  . 

Such  a  process  y( •)  will  be  called  an  observation  semi  martingale. 

bet  f(tj  x(s),  s  <  t)  be  a  non-anticipating  functional  of  x(  •)  such  that 
Bf(t;  x(s),  s  <  t)  <  “,  V  t  8  10, T)  ,  and  define  Fy  s  -  o{y(s)  |0  <  s  <  t}  and 
F*'y  :  -  c{x<s),  y(s)  |0  <  s<  t}. 

Theorem  1.1  (Kalllanpur,  Striebel  (9J ) .  Let 

-  exp  (-  h(s,x(s))dw(s)  -  j  /J  h2(s,x(s))ds) 

Then  (i)  Pn  is  a  probability  measure,  and  P  and  PQ  are  mutually  absolutely 
continuous. 
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(li)  VdP~  1  t'Y}  "  expt^O  h(8,X(8))dy{»)  -  l  fl  h2(s,x<s>)dsl. 


(ill)  On  (Q.Pg),  y ( • )  is  a  Brownian  motion  independent  of  x(>). 

(iv)  x(«)  has  the  same  Law  on  <fl,PQ)  as  on  (ft, P). 

(v)  E {f ( t>x( s) ,  s  <  t)  |  F*} 

E0(f()x(s),s<t) 

-  - — - - -  .  (1.2> 

For  a  nice  treatment  of  this  theorem,  see  Wong  [21].  It  is  the  principal  theoretical 
tool  for  our  work  in  filtering,  for  it  explicitly  characterises  the  optimal  filter  as  a 
functional  integral  and  it  establishes  that  y( •)  is  mutually  absolutely  continuous  with 
Brownian  motion. 

Finally,  we  remark  that  we  restrict  ourselves  here  to  scalar  processes  only  in  the 
interests  of  notational  simplicity.  The  techniques  to  be  discussed  extend  easily  to  the 
vector  case. 


2.  Multiple  Integrals 
2. 1  Definitions 

The  concept  of  a  multiple  Wiener  integral  derives  ultimately  from  Wiener’s  work  on 
’homogeneous  chaos*  decompositions  of  functionals  of  Brownian  motion*  however,  the  modern 
definition  and  theory  are  due  to  Ito  [7] .  Here  we  will  define  multiple  integrals  by 
iteration  of  stochastic  integration.  While  this  differs  from  Ito’s  construction,  it  leads, 
as  Ito  (7)  notes,  to  the  same  result  modulo  a  multiplicative  constant.  The  iterative 
definition  is  convenient  for  our  calculations. 

Imt  (b(t),F^)  be  a  standard  Wiener  process  with  its  associated  family  of 
cr-algebras  F^  *  afb(s)  :  a  <  t).  Recall  that,  for  a  Jointly  measurable,  F^-adapted 
process  $(t,w)  such  that  E  $2(s)ds  <  <»,  the  Ito  integral  $(s)db(s)  has  the 
properties 


♦ 

!•* 


,*r 


E  f*  $(s)  db(s)  -  0  t  <  T 


E(/jj  $(s)db(s))2  -  E  /£  $2(s)ds  t  <  T 


L2((0,T]r)  -  {f  e  L2([0,T]r)  |  f  is  symmetric} 

This  will  be  the  set  of  integrands  for  the  rth  order  integral.  If  f  e  t2l [0,T]r) , 

2  r-1 

f(o, ...)  e  l  (  [0,T]  ),  will  denote  the  section  of  f  at  a. 

Definition  2. t  Let  f  e  L2([0,T]r)  t  <  T.  l£(f)  is  defined  recursively  by  (L2([0,T]°) 


It(f)  «*  f  for  r  ■  0 


if  (f  >  -  f*  Ir_1(f( . .  )db(  s) 

t  ■'OS 


I*<f)  is  the  rth  order  multiple  Integral  of  f  with  respect  to  b( •)  up  to  time 
t.  Alternately  stated. 


xt«f>  “  Jo  Jo^'-Jo"’1  f(S1 . Sr)db(Sr>...db(Sl) 


To  insure  that  the  rlgjit-hand  side  of  (2.3)  is  well  defined  it  suffices  to  show  that 

—  I  2 

Ir-1  (f  ( s , . . . )  )  has  a  jointly  measurable  version  with  bounded  L  (SI  x  (0,T1,  P  x  \)  norm, 
s 

This  may  be  done  by  proving  recursively,  along  with  the  definition,  that 
EI^(f)I^(g)  -  py  (f,g) 

(2.4) 

-  Jo  /o1  ***/or_1  f<V-Sr)9(al . sr)dsr...dSl 

for  all  f,g  e  L2((0,T]r).  This  is. a  consequence  of  (2.2).  Then,  if  fn  is  a  sequence  of 
symmetrizations  of  separable  functions,  such  that  fn  *  f  in  L2-norm,  I* (f n(s, . . . ) )  is 
jointly  measurable  for  all  n  and 


...  •  v 


ft  [I*_1(fn(s,...))  -  I^Vffs ,...))I2ds  -  0  . 


Thus  we  can  find  a  jointly  measurable  version  of  Ig  <f(s, ...)). 

It  is  important  to  note  that  multiple  integrals  have  zero  man  and  that  integrals  of 
different  orders  are  orthogonal!  that  is,  for  f  e  L2((0,T]r),  g  8  L2( [0,Tl**),  q  /  t,  t. 


El‘(f>  -  0 

(2.5) 

Efl*(f)l^(g)!  -  0  . 

These  follow  from  repeated  application  of  (2.1)  and  (2e2># 

Remark.  The  requirement  of  symmetry  for  the  integrands  is  not  necessary,  since  integration 
is  carried  out  only  over  the  set  where  s^  >  >•••>  s^.  However  this  convention  Is 

convenient  in  formulating  the  multiplication  formula  in  section  2*2. 

The  following  technical  lemma,  a  Fubini  result  on  interchanging  db  and  ds 
integrations,  is  needed  later. 


Lemma  2.1  Let  f  e  l2([0,T]r).  For  t  <T 


J*  I^_,(f(s,...)ds  *  /J.../Qr'2  /*  f(u,s1,...sr_1)dudb(sr_1)...db(s1)  .  (2.6) 


Proof  Define 


g«.(s1(...s_  .)  “  f  ( u,  s  ...  *s  ,)du.  The  r.h.s.  of  (2.6)  la  1*  (g  )• 

t  i  r- 1  '8  i  r- 1  v  t 


To  prove  the  lemma,  simply  verify  that 


if/*  Ir'1(f(s,...))ds  -  xrVii*  “  0 
os  r  r 


by  using  the  basic  properties  (2.5)  of  the  multiple  stochastic  integral. 


For  filtering  applications,  we  must  also  define  multiple  integrals 


Jg  •••f0r”1  )dy(8r)  •••(dy(s1) 


with  respect  to  observation  semi-martingales 


4 

4 

* 

* 


y(t)  “  /!(  h(x  )ds  +  w(t) 


(the  assumptions  of  section  1  are  assumed  to  be  in  force).  Such  integrals  are  known  and 
have  been  studied  in  the  context  of  semi-martingale  theory.  However,  the  special  structure 
of  (2.8)  allows  a  simple  definition  which  we  present  here.  This  takes  advantage  of  the 
absolute  continuity  of  the  y(  •)  process  with  respect  to  Brownian  motion;  as  stated  above, 
if  (ft,  F,  P)  is  the  underlying  probability  space,  there  exists  a  probability  measure  PQ 
such  that  P q  <<  P,  P  <<  P0,  and  y(  •)  is  Brownian  on  (ft,  F,  Pq).  Therefore,  for 
f  e  L2([0,T]r),  we  define  (2.6)  as  the  random  variable,  which  on  (ft,  F,  P^)  equals  the 
Brownian  motion  integral  defined  above.  We  call  this  integral  I^(f)  without  reference  to 
measure  or  process,  which  should  always  be  clear  from  context. 

The  iterative  property  of  I*(f)  remains  true  for  dy  integrals;  that  is. 


It<f)  ”  ^0  It-1<f<8....))dy(s) 


where  the  integral  in  (2.9)  la  defined  with  respect  to  the  semi-martingale  y(  •)  in  the 
usual  sense  (see  Liptser  and  Shiryayev  (11]).  However,  neither  the  expression  (2.4)  nor 
the  orthogonality  of  different  orders,  (2.5),  now  holds.  Instead,  we  can  prove  the 
following  lemma,  which  is  useful  in  section  3.2.  (In  this  discussion,  we  abbreviate 
h(s,x(s)>  by  h(s)4 

I,emma  2.2  Suppose  Ef h2(s)ds)r  <  ®.  Then  for  k  <  r  and  f  e  L2((0,T]k) 

(i)  E[I*(f)l2  <  Mklf|22»  <  »  is  independent  of  f 

la 

(ii)  E  I*(f)  ”  /g  1f(s1,...,sk)  E(h(s1  )  •••hfs^)]  ds^  •  •  •ds1 


Proof 


We  will  actually  prove  by  induction  the  more  general  result:  for  r  >  4  >  k 


e£0,T1 


E[x(ot)..-x(Vl)I^(f)]2  <  «U(V, 


(2.10) 


where  a  ,  e  ''l,  and 

1C 


E[h(o  )•••»>< o^)!^  <f)} 
s  °k 


°k  ,S1  f9k-1 


(2.11) 

f0K  /0  1  *  *  *  /0  K  '  f(ai'***,sk>  Eth(81)*»h(8k)h<ak+1)**h(o)l))d8k.*d31  . 

Lemma  2.2  is  the  case  Z  =  k  for  every  k  <  r.  First  we  demonstrate  (2.10)  and  (2.11)  for 
r  >  l  >  k  «  1,  using  the  iterative  formula  of  (2.9)  and  the  Independence  of  x( •)  and 
w(«).  Thus 


E(h(ot)--h(a2)  /0’f(8)dy(a)l2  - 

a1  °t  1 

E(h(aJl)**h(o2){|0  f(s)h(s)ds  +  fg  f(s)dw(s))r  <  (2.12) 

(2E  /jj  [h(al)--h(s)]2ds  +  2  E(h(ajl)2*»h(o2)2]  )lf«2  =  a  1  (  a^,  •  -o^)  If  I2  . 


To  derive  the  Inequality  in  (2.12),  the  Cauchy-Schwarz  inequality  is  used  several  times. 
»t  1  c  L^fO.T)1  ')  far  t  <  r  because  E  ( h2(s)ds]r  <  Likewise 

o, 

EChCOjI-.hi^)  /Q  f(s)dyaJ  - 

o  a.  o 

E[h(  al>  **h(  a2)  [fQ  f(s)h(s)ds  +  fg  f(s)dw(s)l)  **  J0  f  ( s)E  (h(  o2>  •  *h(  c^))  ds  . 


(2.13) 


Now  suppose  (2.10)  and  (2.11)  are  true  for  a  fixed  k  and  all  l,  r  >  l  >  It.  Again, 
using  ij^'ff)  «  Ik  (f  ( s,  •  •)  )dy(s) ,  Cauehy-Schwar  z ,  and  induction 


Elhlo^.-Mo^)!^1  <f)]2  < 

2  C*1  /„Vl  E(h(<j4).*h<alt+2)h<s1)l^  (f(s2,..)))2ds2dBl 


+  2  /Qk  Eth(oe)-*h(o)c+2)lk(f(s,..))12ds 


«  12  ^U(8'V2'"'0t)dl  +  2,H,k(V2',,V1,',J 


ai.k+1(<W'V,f' 


By  induction,  az  fc+1  e  L1  ( (0,T]  *~k_1 ).  Thus  (2.10)  is  true  for  k  +  1.  That  (2.10)  holds 
for  k  also  implies 


E  Ik(f(a,  ••)  )ds  < 


Thus ,  from  (2.1), 


E  f*T  Ik(f  (s,  •  •)  )dw(  s)  »  0,  for  t  K  T  . 
'0  s 


With  the  aid  of  this  equality  we  can  prove  that  (2.11)  also  is  true  for  k+1.  This 

completes  the  induction  step.  Induction  stops  at  k  m  c  since  we  have  required 

,T  2  r 

r  >  *  >  k  in  order  to  apply  E( JQ  h  (s)ds)  <  «. 
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2.3  The  multiplication  formula. 

As  above  let  (b(t),F^)  denote  a  standard  Brownian  motion.  If  ^(b(s),  s  <  t)  is  a 
functional  of  b(«)  up  to  time  t,  we  want  to  consider  expansions  of  the  form 


*  -  l  • 

r-0 

(If  e  L2(n,F^,P)  such  a  representation  exists,  uniquely,  and  the  series  converges  to 
i|>  in  mean-square;  see  Ito  {7}  or  Hida  (5).)  Rules  prescribing  how  this  representation 
changes  as  various  operations  are  performed  on  must  be  available  if  multiple  integral 
expansions  are  to  be  of  use  in  applications.  In  this  section,  we  address  the  simplest 
problem  in  this  direction.  If  f  S  L2([0,T]r)  and  g  6  L2((0,T]'1),  what,  if  any,  are  the 
kernels,  such  that 

I*<f)Iq(g)  -  l  I^U1)  (t  <  T)  7 
c  c  i-0  * 


To  express  the  answer,  we  first  introduce  the  following  notation. 
Definition  2.2 


(1) 


(ii) 


Pr  £  projection  of  L2([0,T]r)  onto  i,2([0,Tjr)j 
<Ptl»<V..,or>  -77^  h'%(,,'**'0x(r)) 


r 

where  Sr  -  permutation  group  on  r  letters. 
For  integers  r,q,k,  0  <  k  <  min(r,q),  and 
f  e  L2((0,T]r),  g  e  L2((0,T]q) 


(f  ■k(t)g)(£r1,..,cr+q_2k> 


HiTT  /o  "Jo  f(V'VV*'°r-k>  g(VV  °r-k+1 '  *  *'  °r+q-2k  ><l8k  *  *  ***81 


(iii)  f  g  =  P^^tf  »k(t)g] 
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<iv)  f  9  g  -  f  ®0<t)g  »  P  (f(  ot ,  •*,qr  )g(  or+t ,  o  )] 


9^(t)  is  the  operation  by  which  new  kernels  are  created  from  old)  indeed. 


f  ®k<t)g  s  L2([0,T]r)  xl2((0,Tlq)  L2(  tO,T]r+q“2k) 


as  the  following  lemma  demonstrates 


Leimpa  2.3  For  every  t  <  T 


In  fact. 


f  »k(t)g  c  L2(t0,»lr+q"2k) 


If  CJt)gl2  <  c  ,  Ifl2  Igl2 
k 


where  cf  k  is  independent  of  f  and  g. 

Proof  It  suffices  to  prove  the  lemma  for  a,  instead  of  9,  since  pr+q_2k  18  8  bounded 
operator.  Let  do  “  doi **dot+q-2k'  d®  ”  d81**d8k*  We  then  have,  using  the  Cauchy-Schwarz 
inequality 

“  'k't’g'  *  / 


,r+q-2k  (kl)‘  [0,T] 


<  Ifl2  Igl2 

(kl) 


To  understand  the  meaninq  of  S^(t),  it  is  useful  to  think  of  the  functions  f  and 
g  as  tensors,  which  they  in  fact  are  by  the  isomorphism 

t2<(0.T]r)  -  L2(tO,T]>  a*’*  L2([0,TJ)  (r-fold)  . 

Then  f  9^1  t)g  may  be  viewed  as  a  tensor  contraction  since  it  ‘sums',  that  is, 
integrates,  f  and  g  along  the  first  k  indices.  Thus  f  9fc(t)g  is  simply  a 
symmetrised,  k-fold,  tensor  contraction.  It  is  in  this  definition  that  the  symmetry  of 
f  and  g  is  used)  otherwise  a(t)  would  have  a  more  complicated  definition.  For 
notational  convenience,  we  shall  often  write  9^  instead  of  9^(t),  in  which  case  the 
(t)  is  to  be  assumed.  When  the  time  parameter  is  important  or  different  than  t,  it  will 
always  be  given. 


I 


.  -  .. 


*-  -ri 


W«  can  now  state  the  result 

Theorem  2.1  Let  f  et2((0,T!r),  g  e  L2(  [0,T] q> .  Then 


min(r,q)  ,  . 

I*(f)  Iq(9)  -  I  I^q  ((rr2)c2,')fk  a  (tig  >  . 

k«0 


(2.14) 


Remarks  1.  (2.14)  shall  be  referred  to  as  the  multiplication  formula.  Rida  has 

independently  derived  this  result  as  an  application  of  his  theory  of  generalised  Brownian 
functionals  (personal  communication  of  T.  Hi da;  for  generalised  Brownian  functional  theory, 
see  Hida  (41).  Our  proof  is  elementary,  using  only  Ito's  differentiation  rule,  for 
similar  theory,  see  also  Meyer  [15].  Versions  of  this  formula  are  also  known  in 
mathematical  quantum  field  theory  (Reed,  Simon  [19]).  See  Mitter  and  Ocone  [17]  for 
further  comments. 

2.  The  multiplication  formula  generalises  a  Hermlte  polynomial  lndentity.  The  nth  order 
Hermite  polynomial  of  a  single  variable  is 


Let 


.  ,  .  (-1)n  2,„  d"  2 

h_(x)  •  — I- —  e  x  /2  -  -  x  /2  . 

IB  rJ 1  B 

dxn 


Gr  -  (I*(f)  |  f  8  L2!  [0,t)r» 


and  let  (4  }"  be  an  orthonormal  basis  of  L2([0,t]).  Then,  (Ito  [7],  Kallianpur  [8]) 

n  n“  1 


_  ■«  t  |  P,  p  -  r  , 

G  -  Sp[  n  h  (/„  ♦.  (s)db(s)  )  .  V  .  ,} 

r  P^O  ■’  I  Ji  are  pairwise  unequal 


where  Sp  denotes  the  closure  in  I.  (»)  of  the  linear  span.  One  then  sees  that  (2.14) 
generalises  the  identity,  ([12]), 


min(r,q) 


hr(x)h(xi-  i  K><q) 


y2 


k-0 


r+q-2k 


(x) 


(2.15) 
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V  r,q  >  0.  There  is  a  discrepancy  between  (2.15)  and  (2.14)  in  the  (actors  multiplying  the 
expansion  terns,  but  this  is  due  to  the  different  normalizations  Involved  in  the 
definitions  of  Zr  and  0.  The  relationship  between  (2.14)  and  (2.15)  nay  be  seen 

clearly  in  Hida's  work,  but  we  shall  not  pursue  the  natter  further  here. 

We  will  show  how  to  prove  theorem  2.1  using  ito's  rule  and  Induction.  For  this 
purpose,  we  need  certain  facts  and  identities  concerning  0,  and  these  are  collected  in 
the  next  leans.  The  notation  f (S^, .. ,3^, ...)  indicates  the  section  of  f  in  which  the 
first  k  variables  are  fixed  at  s  s^  respectively. 

Lemma  2.4 

(1)  (<«,.••)  Ok(c1)g («,.••)  <<V**'°rtq-2k-1)  e  L2(  £0,T)  r4'q“21c) 


(ii)  f  Ok(t)g  -  f  Ok(o)g  ♦  /*  f(st-*)0K_1(s)g(s,..)ds 
(ill)  For  k  >  1,  (f  OjjItlgKSj.'ha^^) 

”  f<  v>  Vt)9  +  f  vt)g(0i '"n 


(02'**'°rtq-2k) 


(2.16) 


(2.17) 


(iv)  (f  0(t)g)(c,,**,c_  )  -  (2.18) 

i  rtq 

(•^  ((c^**)  0(t)g  ♦  f  0<t)q(a1.**)l  to2' *,'°r+q-2k)  * 

Proof  i)  follow  from  calculations  similar  to  those  in  leans  2.3.  The  details  will  not  be 
presented. 

ii)  By  direct  calculation  and  definition,  using  the  symmetry  of  f  and  g  extensively. 

f  Ok<t)g 


-  P 

r+q-2k 

•  P 

r+q-2k 

m  D 

r+q-2k 


(-  ••  J* 

lkl  J0  J0 


|f*  f’1... 

lJo  Jo 

c—  fa  •• 

lki  J0  J0 


ds1*«dsk  f(s),»»,sk,»*)g(s1. 


/,•  f(a1.**,ak,*»)g(s1,**,sk. 


•)dsk  ••dS1) 
• )  dsk  •  *ds  ^  1 
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I 


+  Pr+q-2k{(1^7TT  £  /o1"  C'1f<V**'V‘*>9<V*,'V">dV,,ui1 

-  f  »k(o)g  +  J*  ds  f(s,**)  9)(_1  (s)g(s,  ••) 


(ili)  and  (lv).  The  proof*  of  (ill)  and  (iv)  are  sisdlar,  (lv)  being  just  a  special  case 
of  (iii).  He  ahall  only  present  (iv),  as  it  is  simpler.  Note  first  that,  by  definition, 

—jlf («,»••>  «  <t)gl  <02'**'°r+q> 

”  r+q  ( r+q- 1 )  I  J-  f  1  V  °*(2 ) '  * ‘°»(r ) )  9<  °i(r+1 ) '  * *'  °ii(r+q) )  (2.19) 


where  w  e  S  .  is  interpreted  as  a  permutation  of  {2,  •  *,r+<i).  Now  using  the  syaswtry 
r+q- 1 

of  f,  (2.19)  may  be  written  as: 

,  r 

J-  Ua*i2)'”'e«i)'  ®1'®t(J*1)'*,'*w(r))  * 


l-1  ¥E£Wi 


(2.20) 


9<0t(r+1)',*'0i.(r+q)>  * 


Using  the  expression  analogous  to  (2.20)  for  ■jJL  f  a  (t)g(  o1  ••), 


f(o,  ,••)  9(t)g  +  f  «(t)  g(o  ,  ••)}  (a **,0  ) 

r+q  1  r+q  i  i.  z+q 


mT^knlJ-  ml  °s(2) '  **' ai' **' °*(r)  >  *  °*(r+1) '  *  *' °»( r+q)  ^ 


1-1  weSr+g-l 


+  A  J  f*°w(2)'**'°w(r+1)*  X  ,<<,w(r+2)'**'01,**'°»(r+q))1 


j-1  ires 


Tr^qTT  .1  f(0x(1)''*)9‘0ir(r+1>'-,'°r(r+q)) 


■  f  9(  t)g(  a, ,  •  • ,  o  ) 
1  r+q 


"  v<V  - .  *•* 
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This  la  tha  desired  raauXt 


Proof  of  theorem  2. 1.  Ha  use  Ito’a  differentiation  formula  and  the  preceding  lemmas  to 
implement  an  induction  argument  that  proceeds  in  two  steps i 

(a)  Show  (by  induction)  that  (2. 14)  holds  for  orders  r  “  n,  q  ■  1»  Vn. 

(b)  Assuming  (2.14)  for  (r-1,q),  (r,q-1)  and  (r-1,q-1),  show  that  it  holds  for  (r,q). 
(a)  and  (b)  than  provide  a  consistent  scheme  of  induction  for  proving  theorem  2.1  for  all 
orders. 

Step  (a)  By  Ito's  differentiation  rule 


j||  f(s)db(s)  fjj  g(s)db(s)  -  /Jj  /0’  [f(s1)g(s2)  +  f(s2)g(s1  )]db(s1 ) 
+  f(  s)  g(  s)  ds  . 


This  proves  the  case  r  *■  q  *  1. 

Suppose  that  the  theorem  is  true  for  (r,q)  -  (n-1,1)  and  let  f  t  L2(I0,T]n), 
g  c  L2((0,T] ).  Applying  Ito’s  differentiation  rule  again. 


l£(f)l’(g)  -  J*  g(s)  lj|<f)db(8>  +  /*  i“_1  (f(s,..))i’<g)db(s) 


(2.21) 


By  induction, 

In-,(f(s,.*))ll(g)  -  l"(n(f(s,..)  0  g]>  +  l"”2(f(s,..)  O  (s)g>  . 

(  8  8  a  i 

Lemma  2.5(1)  and  lemma  2.1  justify  interchanging  integrations  in  the  last  term  of  (2.21): 


ft  l""\g(a)f(s,.0)ds 


Thus,  by  substitution  in  (2.20) 


l£(f>  l’(g>  -  /g  {l^( g(  s)f  {  ••) )  +  I™(n(f(s,  ••)  0  g])}db(s) 

+  /g  ei(8)g)db(a) 

+  l"  1 1  g(u)f(u,«1 ,  )du) 


-  l£  {g(c1)f(<»2.**«on)  +  og< «1 )]  (Oj»  “.on) } 

+  l"  ©1<o1)g<o2,  •**«„>  +  /q  g(s)f(s,  Oj,  •♦,an_1)ds}  « 


By  lemma  2.5  (ill)  and  (iv)  this  becomes 

l"+1  ( (n+1)f  0  g>  +  l£-1  (f  31(t)g)  , 

which  completes  the  Induction  step  of  (a). 

Step  b  Without  loss  of  generality  assume  that  q  <  r.  The  Induction  hypothesis  is  that 
theorem  2.4  is  true  for  (r-1,q),  (r,q-1),  and  (r-1,q-1).  Apply  Ito's  differentiation 


I*<f)I*(g>  -  Jg  l^(g)I^'1(f(s,*.))db(s) 

+  I?"’ <9)ir<f )db(s> 

'U  B  8 

♦  /g  1(f<s,*.))I^"1(g(s,.»))ds  . 

Next,  use  the  Induction  hypothesis  to  expand  the  integrands  in  (2.22),  then  interchange 
da  and  db(s)  integrations  where  necessary,  and  collect  like  order  terms.  The  result  is, 
for  q  <  r 
I*(f)l’<g) 

”  {t(r^11)lf(s1,»»)  0  gl  +  ir+^"1)(f  0  g(s1,**)ll(s2,*«,s  )) 


>  (^*2k)  tf  0k(.1)9(-1,-.)1(-2,...sr+q.2k) 


-,v 


+  <r^]c21')  J*  f(u,  ••)  e)t_1(a)g(u,  ••)du} 


♦  I 


^  q{{f<«1,  ••)  e  (a^g]  («2<  •*,«  )  +  /*  f(u»  ••)  O  ^(u)g(u,  »»)du) 


To  complete  the  proof,  we  need  only  apply  the  identities  of  lemma  2.4  (Hi)  and  ( iv)  to  the 
kernels  of  this  last  expression.  For  example,  the  kernel  of  I^+q-2k,  1  <  k  <  q-1 ,  equals 

(r+r-k2k)  t~Ek<f(V>  V-1>q)<*2'”>  ‘V> 


+  (/*  f(u,  ••)  °K<u)g<u,  «*)du)  (s2,  ••)) 


l(f  0K  («)lq!(»1,,,,«  r+q-2k>  +  (  *<«#••)  Ok_,  (u)g(u,  •  »)du)  (Sj,**) 


-  Ckf,f(f  Vt,q)(V'Vq-2k>5  • 


This  is  the  kernel  given  in  (2.14).  The  kernels  of  I^*1*  2*/  k  »  0  and  k  -  q  are 
treated  similarly.  This  completes  the  proof. 


3.  Multiple  Integral  Expansions  in  Filtering  ‘Rieory. 

This  section  explores  the  use  of  multiple  integral  expansions  for  optimal  and 
suboptlmal  filtering.  The  estimation  problem  considered  is  the  general  problem  stated  in 
the  introduction,  and  the  notations  and  assumptions  established  there  shall  remain  in 
force.  For  additional  notational  convenience,  we  let  f(t)  :  -  f(ti  xg ,  s  <  t),  h(s) 

:  -  h(s,  x(s))  and  ffc  s  -  E{f(t)  |  F^}. 

3. 1.  Expansion  of  the  optimal  filter 

A 

In  theorem  3.1  below  we  derive  an  expression  for  ft  as  a  ratio  of  two  multiple 
integral  expansions  in  which  the  process  of  integration  is  y(t),  the  observation  semi- 
martingale,  and  the  integrands  are  deterministic  functionals  computable  from  the 
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C~ - s - - - 


(unconditioned)  dietribution  of  the  signal  process.  First  we  state  some  preliminary 
definitions  and  a  lemma. 

Let 

Lfc  t  -  exp{ J*  h(s)dy(s)  -  Vi  Jg  h2(s)ds]  . 

sp 

Lt  is  the  important  process  in  this  calculation.  Observe  that  ty  -  and 

0 

(Lt#  F*  v  F*)  is  a  martingale  on  (0,  F,  Pfl),  (F*«  “  o{x(s);  s  e  R+)>.  h  conditioning 
argument  then  shows  that  the  Kallianpur-Striebel  formula,  (1.2),  can  be  expressed  as 


*0{f(t)Lt  I  I») 
f  _  - - - 

Vlt  1 

The  following  process,  based  on  1^. ,  will  also  appear t 


(3.1) 


,  <  r) 


H  h« 


1>....h(sr)dy( sr>. 


•Ws,) 


Note  that  L^r*  la  not  a  multiple  Integral  of  the  type  defined  in  §2  since  the  integrand 
is  not  deterministic.  L^r*  may  be  properly  defined  by  noticing  that 

.  (r)  ,t  .  (r-1 ) .  ,  . 

-  /0  h(s)Ls  dy{  s)  . 

Iterative  use  of  the  stochastic  Ito  Integral  then  specifies  L^r*  for  any  order  r.  This 
is  especially  easy  to  carry  out  on  (ft,  F,  PQ ) ,  on  which  y(t)  is  a  Brownian  notion 
Independent  of  the  signal  (see  theorem  1.1). 

The  following  stochastic  Pubinl  theorem  for  interchanging  conditional  expectation  and 
stochastic  Integration  is  needed;  lt  is  a  direct  consequence  of  theorem  S.14  in  Llptser  and 
Shlryayev  (11). 

Leiwsa  3. 1  Let  *(  s)  be  a  Fxy  adapted  process  such  that 
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>*Ar  -u  v 


VG  *2(s)da]  <  «  . 

Than  E{^  *(s)dy (a)  |  F^>  -  Jg  Sg(^B)  |  F*}dy(a). 

Finally,  it  ia  convenient  to  introduce  the  functiona 

i  (t,  a, . .. ,a  )  i  -  E{f (t)h(a  >...h(a  >}  n  >  0 
n  l  n  in 

k  (a.,. ..,a  )  t  -  E {h( a  ) . • .h( a  >}  n  >  1 

n  i  n  i  n 

kQ  :  -  1  . 

Theorem  3.1 

i)  (Partial  expanaion)  If  *(/g  h2(o)do)r  <  ”,  and  E[f2(t)(/p  h2(o)do)rl  <  •», 

then 


I  I^n)<tB<t>>  +  Efl{f(t>  L**'  |  F 

n-0  _  _  _  _ 


tr>  i  rn 


(3.2) 


i  i;n,(kn>  ♦  *0{f(t)  L*r)  i  f*> 

n-0 


11)  (Full  expanaion)  If  E(exp  h2(a)da)  <  ”,  and  E(f^(t)  exp  h4(B)da]  <  ” 


•t  .2, 


I  i[n)un(t>) 


2  _  n-0 

ft 


(3.3) 


i  ir,(kn> 

n-0 


and  the  expanalona  converge  in  L1  (P) . 

Proof i 

Part  i)  By  applying  Ito'a  differentiation  rule  to  , 


dLg  -  h(a)Lsdu(a) 


mo  that 


-  1  +  h(*)L  dy(a)  . 


(3.4) 
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Iterating  (4.4) ,  we  find  that  for  any  r 


I*t  “  1  +  /g  h(s)dy(s)  ■*•  fa  h(.1  )h(»2)dy(s2)dy(s1 ) 


+  • • ••  +  L 


Mow  substitute  this  expansion  into  the  Kailisnpur-Striebel  formula  3.1  for  ft.  The 
denominator,  for  example,  becomes 


WFt}  “  1  +  l  Eot/o*’*^on1h<,1)*“h<sn)***dy<81)IFt} 

n-1 

+  *o{Ltr>1Ft}  * 


The  hypothesis  E[ h2(s)dslr  <  »  of  (i)  allows  lesssa  3.1  to  be  applied  to  the  terms  of 
(3.5),  with  the  result. 


*0[i-t|Fyl  -  1  +  l  /g*««/on-1E(|{h(s1)  •••h(sn)}dy(sn)  •••dyts^  +  E{L^r)|Fy} 


(r)lty. 


Since  the  distribution  of  the  signal  process  is  invariant  under  the  change  of  measures 


from  P  to  P„ 


Therefore 


E  {h( s. )  •  •  *h(  s  )}  »  E{h(s  )  •  •  »h(s  )} 
u  i  n  in 


»  h  (s  ,...,s  ) 
n  1  n 


“0{tt,F^}  *  l  X[B>«V  +  E^r>IFt)  * 


h  similar  calculation  yields 


EQ{f(t)LtIFy}  -  l  l£n)U  (t))  +  E {f ( t )L^r)  IFy } 

n«0 


Substitution  of  these  expressions  into  the  Kallianpur-Striebel  formula  then  proves  (3.2). 
Pert  ii).  Formally,  the  proof  of  the  full  expansion  follows  by  setting  r  •  >.  To  prove 
it  rigorously,  we  first  show  that  E  exptj^  h2(s)dsl  <  •»  implies 


m.s.(P0)  lim(1  +  1  /£•••/, 

N*<»  n-1 


n-1 


h(s  )*»»h(s  )dy( s  )***dy(s,)J 
1  n  n  I 


(3.6) 


Denote  the  finite  series  on  the  right  hand  side  of  (3.6)  by  a”.  Then 


WAt>2 "  Eol^o"*Ch(8i)"*h<sN+i)Ls  dy(8m-i)‘*,dy(8i)]‘ 

N+ 1 


By  employing  the  standard  computational  rules  (2.1),  (2.2)  for  stochastic  integrals,  this 
last  expression  equals 


s 

/S—Cv 


v, 


provided  that  it  is  finite.  However, 


E0(h2(s1>...h2<sN)L2 


N+1 


EoO''C(s1)**’h^t8N)  expl-/0N+1h2(s)ds)  E(j(exp[2  /QN+1h(  s)dy( s>]  IF*  ]  }  .  (3.7) 


Now  on  (Q,P  ),  x( •)  and  y(  •)  are  independent  and  y(  •)  is  Brownian,  and  hence,  given 

9 

(x(s),  s  <  Sj),  f0  h(s)dy(s)  is  a  Gaussian  random  variable  with  mean  0  and  variance 

rSN+1.2,  _ 

h  (s)ds.  Thus 

S  S 

EQ (exp  2  /gN+1h(  s)dy(s)  IF*  )  -  exp  2  /g^h^slds  .  (3.8) 
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Therefore,  using  (3.8)  in  (3.7) 


(3.7)  -  £0{h2(s1)*»*h2(stJ+1)  exp(/0l,+1h2(s)ds3  } 


e0{hZ(s1)...h;4(stI+1)  l  1~V(o1>***h8(oj>a«J—d«|)  . 


As  a  result 


/o—/oR+lBolh2<8i)— h2(8Hti>L^  ,ldvr*‘i 

NtI 

“  I  /o*,*/oj"lEoCh2<8i)“*h2'8j>ld8j‘“d8i 

1 


J+i  JT  h2(s1)...h2(sj)dsj...ds1} 


(3.9) 


Since  E  exp[  J^h2  (s)ds]  <  «,  (3.9)  tends  to  0  as  N  ♦  proving  that  Lfc  - 

m.s  (P„)  lim  A N  for  all  t  <  T,  as  desired.  Lemma  4.1  can  now  be  invoked  for  every 


order  n,  so  that 


V\|^>  ”  E0{ms*  iim  At^t^ 

N 

“  m.s.  lim  e0{a"|f£} 
N+*» 


.s.(P  )lim[1  +  l  l"(kn)) 
N-h»  n-1 


A  similar  proof  expands  EQ  (f ( t )Lt | F^}  in  the  series 


V«  +  I'l'V  • 

n-1 


Finally,  to  derive  the  L1 (P)  convergence,  note  that 


E  ]  2 
0  dPfl 


EqL2  ■  Etexp  /Jh2(s)ds]  < 
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(ElVVF£5  -  (1  +  [  xt<V') 


<*0  ^,2  E0  -  <'  +  JIt<kn,))2 


1  N 

Thus,  from  (3.6),  E-tVlF*]  -  (L(P))  lim  [1+7  l"(k  )]  as  claimed.  This  completes 
D  t  t  u.  t  n 


N+®  rv»  1 


the  proof  of  theorem  3.1. 


Let  PU,t|F*>  “  E[1(x(t))| F%)  denote  the  conditional  distribution  of  x(t)  given 
tat 

the  observation  up  to  time  t. 


Corollary  3. 1  If  K[exp  /Th2(s)ds]  <  • 


p(  A,t|F*) 


El  ( x( t) )  +  l  I*)(E1  (x(t))h(s  )  •••h(s  )) 
A  "  t  A  i  n 


1  +  I  l"  (Eh(s^)  •••h( sn> ) 
n-1 


A  related  formula  is  also  of  interest.  If  x(t>  has  a  density  q(x,t),  x(t)  has  a 
conditional  density  given  by 


E-OU  t)  |F^,x(t)-x)q(s,t) 
p( x,t |F^>  *  — - S - - -  . 

^  m  r*  f  A.  \  I  rJ  t 


E0tL(t)|F^) 


Using  the  same  techniques  as  above,  we  can  easily  derive 


E0(I.(t)|F^,x(t)-x}q(x,t)  -  (1  +  J  l"(E[h(s1)..*h(sn)|x(t)-x]  )) 

n^O 

*  q(x,t) 


(3.10) 


for  the  numerator  of  p(x,t|F^).  (3.10)  is  often  called  the  unnormalized  conditional 

density. 

Remark:  These  results  all  have  an  obvious  generalization  to  the  multidimenslnal  case. 

A 

The  Bayes  formula  (4.1)  for  f t  is  properly  viewed  as  the  ratio  of  two  conditioned 
functional  Integrals,  in  which  the  dependencies  between  x(»)  and  y( •)  are  linked  in  the 
I»t  term.  The  expansions  of  theorem  4.1  in  effect  calculate  these  functional  Integrals  by 
expanding  1^.  The  x(*)  and  y( •)  interactions  are  then  separated  in  the  sense  that  the 


"...  1  i  ’ — .■  w.  v.'1  ‘  i" •  .-a 


calculation  of  the  filter  ia  decomposed  into  two  parts?  first,  computation,  off-line  and 

prior  to  filtering,  of  the  kernels  l  and  k_,  and,  second,  stochastic  integration  of 

n  " 

these  kernels  against  the  observations.  Of  course,  in  actual  practice  one  can  only  compute 
a  finite  number  of  terms.  In  fact,  if  the  kernels  are  separable  or  are  approximated  by 
separable  versions,  a  truncated  expansion  may  be  realised  in  a  finite  dimensinal  and 
recursive  manner,  because  a  stochastic  differential  system  can  be  constructed  to  realise 
any  multiple  integral  with  a  separable  kernel.  However,  caution  must  be  excerclsed  in 
approximating  the  optimal  filter  by  truncations  in  (4.3),  because  truncation  of  the  series 
in  the  denominator  can  be  a  source  of  severe  instability.  Although  H{L(t)lF^}  >  0  a.s.,  a 
truncation  approximation  may  pass  through  0  and  ro  lead  to  s  singularity  of  the  filter. 
Thus  an  Independent  estimate  of  the  denominator  is  in  general  required. 

Recently,  attention  has  focused  on  the  unnormalised  conditional  density  and  the 
corresponding  'unnoraalised'  conditional  momenta,  which  are  just  the  numerators  of  the 
Kail lanpur-Str label  formula.  B.  Hong  [22]  has  given  a  class  of  Harkov  signals  for  which 
analytic  expressions  of  (3.10)  are  available.  Again  truncation  of  (4.10)  will  in  general 
yield  functions  that  attain  negative  values,  for  this  reason,  cuimulant  expansions 


.  f  l"(g  (x)) 
,  . >  o  t  ’n 

p(x,t|F*)  -  e 


have  been  studied  as  an  alternate  source  of  approximate  filters  (see  Kterno  [1]).  We  will 
not  pursue  these  issues  further,  but  instead  turn  to  other  theoretical  developments  based 
on  therem  3.1. 


3.2  Beat  rth  order  filters 

finite  sums  of  multiple  integrals  provide  a  natural  class  of  causal  functionals  for 
the  design  of  suboptimal  filters.  He  introduce  the  following  definition: 


Definition  3.1 


ii)  The  best  rth  order  estimate  of  t  ( ■)  at  time  t  la  an  element  ?(t)  6  Y^tt)  such 
that 

B(f(t)  -  ?(t>>2  <  E( f ( t )  -  b(t))2  (3.11) 

for  all  b(t)  e  Yr(t).  The  kernels  of  f(t>,  denoted  by  a„(t),  a1(t),...,ar(t).  are 
called  the  optimal  kernels.  A  process  ?(t)  e  Y^ft),  t  <  T,  satisfying  (3.11)  for  t  <  T 
Is  called  the  best  rth  order  estimate  of  f( •). 

Notice  that  the  best  1st  order  filter  is  simply  the  linear  filter,  and  thus.  In  the  context 

of  multiple  integral  expansions,  best  quadratic  (2nd  order),  cubic,  quartlc,  etc.  filters 

are  the  natural  extensions  beyond  linear  filtering. 

In  this  section  we  characterize  the  set  of  optimal  kernels  as  the  solution  to  a  system 

of  linear  Integral  equations.  The  construction  of  these  equations  and  the  proof  of  their 

validity  utilize  the  expansion  formulae  of  theorem  3.1  and  the  multiplication  formula  for 

multiple  integrals  of  theorem  2.1.  Suppose  for  the  instant  that  the  full  expansion  (3.3) 

holds  for  the  optimal  filter  and  that  f(t)  -  J  I?(an(t))  is  an  element  of  Yr(t),  not 

n-0 

necessarily  the  best.  If  f(t)  is  to  be  a  good  approximation  of  f(t),  we  want 


f(t)  -  f<t) 


l  I^(t.tt)) 

1-0  3 


or 

f(t)  l  I3(k  >  -  l  I3U  (t)) 
j-0  3  3-0  3 


(3.12) 


Now  notice  that  the  left  hand  aide  of  (3.12)  can  be  rewritten  as  a  multiple  Integral 
expansion  by  applying  the  multiplication  formula.  In  fact 


f(t)  l  l3(k  )  -  l  l3(g,(t)> 
j-0  e  3  n-0  3 


v”  •,  *  -  <r;l>  °i<t>  kn 

J  (n,n,i)eC^ 


(3.13) 
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where  Cj  >  “  {(m,n,i)  |m+n-2i  -  J,  1  <  min(m,n),  ■  <  r).  Thue,  one  way  to  pick  an 
approximation  f  would  be  to  chooee  the  kernale  an(t)  so  that  g^(t)  matches  VW 
for  as  many  orders  j  as  possible.  In  fact,  this  Is  a  prescription  for  ths  optimal 
kernels. 

Theorem  3.2.  Assume  E(/^h2(s)ds)2r  <  ••  and  E  f2( t)  (J^h2( s)ds)2r  <  »  . 


r 

Then  a  best  rth  order  estimate  exists.  It  is  given  by  f(t)  “  £  l"(a  (t)>  iff 

n-0  t  B 

gj(t,st . . ^ )  «•  E{f(t>h(s1>...h(sJ)} 

<  •  l^(t,s1,...,s^)) 


(3.14) 


for  0  <  j  <  r. 

Remark.  The  equations  at  (3.14)  comprise  r  +  1  integral  equations  for  the  r  ♦  1 
optimal  kernels  a^(t)  0  <  j  <  r.  This  can  be  seen  from  the  definition  of  gj(t)  and  A 
and  will  be  illustrated  explicitly  in  the  examples  to  be  discussed. 

Before  proving  theorem  3.2,  we  first  establish  scan  preliminary  lemmas.  The  first 
deals  with  existence  of  estimates. 

Lemma  3.2.  If  E( h2(s)ds]r  <  «,  then  the  best  rth  order  estimate  exists  and  is  unique. 

k  2  2 

Proof  From  lemma  2.2  E[It(s)l  «  M^lal  j  for  k  <  r.  Therefore  Yr(t)  is  a  mean- 

L 

square-closed  (Hilbert)  space  of  random  variables.  The  lemma  follows  by  the  projection 
theorem. 

Of  the  next  two  lemmas,  the  first  introduces  the  optimal  estimate  to  compare 

suboptimal  estimates,  and  the  second  verifies  a  technical  identity. 

2  v 

Lemma  3.3.  Let  z ,  v  e  L  (!J/',P).  Then 

E( z  -  f(t))2  <  E( v  -  f( t) ) 2  iff  E(z  -  f(t>)2  <  E(v  -  f(t))2  . 

Proof.  Simply  note 


E(z  -  f(t))2  -  E(z  -  f(t))2  +  2E(  z  -  f(t))(f(t)  -  f  ( t) ) 
+E(f(t)  -  f(t))2 

-  E(z  -  f(t))2  ♦  E(  f  ( t )  -  f(t>)2  . 
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Proof  of  theorem  3.2  Mctuat  of  1< 


3.3  it  tuff too*  to  show  (3.14)  holds  If  and  only  If 


S(f(t)  -  f(t)]2  <  E[c(t)  -  f(t)]2 
for  all  c(t)  e  »r(t).  Sinca 

*lc(t)-?<t)l2  -  E(c(t)-f(t)]2  +  S[?(t)-f(t)]2  +  2E[c(t)-f(t)l  ( f<  t)-f(t)  J 
this  will  occur  If  and  only  If 

E(c(t)-f(t)]  [f(t)-f(t>]  -0  V  c(t)  tytl  . 

Thus,  m  will  demonstrate  (3.20).  Begin  by  noting  that 


dP. 


<V%'Ft’> 


(e, 


yiv-1 


,<W> 


Then 


E[c(t)-f(t)][f(t)-f(t)]  -  El¬ 


iot  t)-f{t))lfU)B0{LtlF*)H^{f(t)LtlF*)] 


VVFt] 


*{K(^-,F^l(c(t)-f<t)>(?<t)«0CttIF^J  -  *0{f(t)LtIF^]} 


E0{(c(t)-f(t))(f(t)E0tl.tIF^)  -  E0tf(t)LtIF^l  )}  . 


Next  note  from  (3.13)  that  g^(t)  depends  on  kernels  kn  of  at  most  order  j  + 


f‘t)»o»tl  t>  ■  lj0  1J(kJ)  * 

-  I  l?(g.(t))  ♦  l  I*(g,(t>>  +  ?(t)E  {L*2r>|F*} 
1-0  3  j-r+1  3 


where  the  g^(t),  r  *  1  <  j  <  3r  are  determined  by  the  multiplication  formula, 
partial  expansions  of  theorem  3.1,  we  then  see  that  the  expression  T(  t)EQ  (I<t I F^ 
E0(f(t)LtIF^}  appearing  In  (3.21)  equals 


(3.20) 


(3.21) 


r.  Thus 


Using  the 


<3.22) 


c  2r 

I  i*<g,<t>  -  t.(t>>  ♦  l  i*(g.(t>  -  t,(t>> 
j-0  3  3  j-r+1  3  3 

♦  ^  *  ?<t)»  {L‘2r)|F^> 

j-2r+1  c  3 

-  K0{f<t)L‘2r)|F£>  . 

Since  y(  •)  is  Brownian  on  (fl,PQ),  multiple  integrals  of  different  orders  are  orthogonal 
on  (fl,PQ),  and  so  if  (3.22)  is  ussd  in  (3.21)  we  find 

(3.21)  -  [c0(t)-a0(t)l  (g0(t)-l0(t)J  + 


+  B0{(c(t)-f(t))f(t)*0(L^2r>|F^J  } 

-  B0{(c(t)-f(t))K0(f(t)L‘2r)|F^n  ,  (3.23) 

The  last  two  terms  of  (3.23)  are  zero  by  lemma  3.3.  Thus,  it  is  clear  that  (3.23),  and 
hence  (3.20),  is  zero  iff 

*  ‘j  0  <  j  <  r  . 

This  completes  the  proof . 

The  technique  of  theorem  4.2  extends  to  other  problems  as  well.  Suppose,  for 
instance,  that  a  filter 


a'(t)  -  a'( t)  ♦  \  I^(a;(t)) 

J-1  3 

of  order  q  is  available!  a'(t)  need  not  be  the  best  qth  order  filter.  Let  r  >  q, 
and,  rather  than  ash  for  the  best  rth  order  filter,  let  us  seek  the  "best  rth  order 
corection*  to  a'(t),  i.e.,  the  mean-squara  minimizing  a(t)  of  the  form 
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r  . 

a(t)  -  »'(t)  +  l  I3(a  At)) 
j-7+1  *  3 

where  a^( t) ,  j  ■  q  +  are  free  to  be  choeen.  Define  the  kernels  9j(t)  as 

before,  but  with  a^(t)  replaced  by  a^(t)  for  0  <  j  <  q. 

Theorea  3.3.  Let  the  hypotheses  of  theorem  4.2  hold.  Then  a(t)  is  the  best  rth  order 
correction  to  a'(t)  if  and  only  if 

g^(t,Sj ,  •••a_j>  »  B{f(t)h<s1)»**h(Sj)>,  q  +  1  <  j  <  r  .  0.24) 

Proof.  As  before,  it  suffices  to  show  that  (3.24)  holds  iff 

E(c(t)-a(t)l  la(t)-2(t)]  -  0 

r  , 

for  all  c(t)  -  a*(t)  +  £  rlc.ltl).  By  the  same  calculations  as  in  theorea  4.2 

j-qtl  *  3 

E[c(t)-a(t)]  (a(t)-f(t)] 


-  B0{(c(t)-a(t)][a(t)E0{Lt|F^}  -  E0{f(t)LtIF^})} 


•  v  l  l  ijt91<t)-t.(t»  ♦  i  ij(q,-t.) 

j-q+1  3  5  j-o  3  3  j-r+1  3  3 

+  l  I^(q.(t))+a(t)En(L'2r)Fl[]  -  E  {f(t)L*2r)  |F*}1} 

J-2r+1  *  3  Otto 

r  t  *1-1 

“  £  4*"V  (cJ-.j)(t,.1,...,sj)IViJJ(t,.1,...,sj)d.j...d.1 


This  equals  aero  iff  g^  -  for  q  ♦  1  <}  <r. 

Remark  Clearly,  an  analogous  result  holds  for  the  case  in  which  an  arbitrary  subset  of 

(a^ is  given  and  the  remainder  are  chosen  as  to  optimise  the  mean-square  filter 

error.  Thus,  if  a,,  j  e  (i,  #•••»!  }  c  (0,1, •••,r}  are  given,  then  the  (a.(t)>, 
j  i  q  —  j 

J  t  (j.,,,,»j  }  are  optimally  chosen  iff  gj  -  l.  for  every  j  e  {0, 1,*««,r)  - 
1  q 

(j,.--jq}. 

As  a  first  example  of  theorem  3.2  let  us  compute  the  kernel  equations  for  the  best 
linear  estimate  f(t)  -  aQ(t)  ♦  /ga^t.sjdyl  a) .  Prom  (3.13), 


31- 


g„(t)  -  «0(t)  +  /g»1(t,o)E(h(o)do 


9|  (t>s)  -  +  j£a1<t,o>E[h(s)h(c>]do  ♦  aQ(t)Kh(  •)  . 

The  kernel  equations  are  then 

a0(t)  +  /jja^t.oJEMoldo  -  Ef(t) 


a0(t)*h{s)  -  aj(t.a)  ♦  ^at(t,a)l(h(«)h{o)]d9  -  Ef(t)h(s)  , 
or,  eliminating  ag( t)  from  the  second  aquation, 

-0<t)  +  /Ja^  (t,o)B|h(  a)]  do  *  Ef(t) 

a1  (t,s)  +  /^(t.oJcovChls)  ,h<  o)]do  -  cov[f(  t)  ,h(s)  ] 


(3.25)  Is,  of  course,  the  well-known  Wiener-Hopf  type  equation  for  optimal  linear 
filtering.  Before  examining  higher  order  examples,  we  will  discuss  the  Kalman  filter. 


3.3  The  Kalman  filter 

Consider  the  filtering  problem  in  which  h(t,x)  -  H(t)x  and  x(t)  is  a  Gauss-Msrkov 
process  arising  as  the  solution  of  the  system 

dx{ t)  -  Kt)x(t)dt  +  G(t) db( t) 

where  Xg  “  constant  or  a  Gaussian  auv.  independent  of  the  Brownian  motion  b( ♦).  The 
celebrated  Kalman-Bucy  theorem  states  that  the  optimal  state  estimator  S(t)  - 
*{x(t)|F^}  satisfies  the  equation 


dx(t)  -  ?(t)x(t)dt  ♦  P(t)HT(t)(dy(t)  -  H( t)x( t)dt] 
'x(O)  -  Xg 


(3.26) 
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where  P(t)  is  the  solution  of  a  deterministic  Riccati  equations*  It  follows  that  x(t) 
is,  in  fact,  a  linear  functional  of  y( •)?  if  4(t,s)  denotes  the  state  transition  matrix 
of  F(t)  -  P(t)HT(t)H(t),  the  solution  to  (3.26)  is 


x(t)  -  *<t,0)xQ  +  t(t ,s)P(s)HT( s )dy( s )  . 


(3.27) 


This  simple,  linear  structure  is  not  an  immediate  consequence  of  the  expansion  formulae  of 
theorem  3.1,  because,  even  in  this  case,  both  numerator  and  denominator  series  will  be 
truly  infinite  sums.  It  is  therefore  of  interest  to  see  how  x(t)  can  be  derived  from  the 
general  expansion.  He  will  show  that  this  can  be  done  using  theorem  3.2  and  moment 
equalities  for  Gaussian  random  variables. 

The  most  common  proof  of  the  Kalman-Bucy  filter  invokes  the  stochastic  differential 
equation  for  the  conditional  moments  (cf.  Fujisaki,  Kallianpur,  Kunita  [2]).  In  this 

a2  'S 

approach,  the  equation  for  x(t)  requires  knowledge  of  x(t) ,  that  for  x(t)  knowledge 

A, 

of  x  (t),  and  so  on,  thus  leading  to  an  infinite,  coupled  set  of  equations.  To  derive 
the  Kalman-Bucy  theorem,  it  must  be  independently  argued  that  the  conditional  distribution 
of  x(t)  given  is  Gaussian.  Because  of  identities  between  different  moments  of 

Gaussian  m.v.'s,  this  allows  the  moment  equations  to  be  truncated  at  n  ■*  2  and  leads  to 
(3.26)  and  (3.27).  By  way  of  contrast,  the  derivation  here  will  not  require  explicitly 
knowing  the  conditional  density.  For  other  methods  of  deriving  the  Kalman-Bucy  filter,  see 
Van  Schuppen  (20). 

In  the  interest  of  computational  simplicity,  we  will  consider  only  the  most  simple 


dx( t)  -  db( t)  x(0)  -  0 


dy(t)  -  x(t)dt  +  dw(t)  y(0)  -  0  , 


(3.28) 


where  b( •)  and  w( *)  are  independent,  standard  Brownian  motions.  The  techniques  work 
also  for  the  general  case. 


Theorem  3.4  x(t)  ”  j!j  a(t,s)dy(s)  where  a(t,s)  satiefiee  the  Wiener-Hopf  equation 


a(t,s)  +  a( t ,  o)  min(s,  o)do 


a  t  >  a  . 


Before  presenting  the  proof  we  must  recall  the  following  aoaent  identities  (Miller  [16], 
Marcus -Willsky  [14]). 

Lemma  3 . 5 v  Let  (z zk]  be  a  jointly  Gaussian  random  vector.  Then 


B(s  1  -  JiU  —I  +  l  eovfz  ,z  JEf  K  i  ] 

j-2  3  fclj  * 


Proof  of  theorem  3.4.  Since  y(  • )  is  continuous  and  Gaussian,  the  set  of  polynomials  in 
2  v 

y(  •)  is  dense  in  L  (8,F*,P),  (Kallianpur  [8]).  Therefore,  it  suffices  to  show  that 
a(t,8)dy(s)  is  the  best  rth  order  estimate  for  every  r,  1  <  r  <  ••  Since 


E <£  b2(s>ds)r  < 


E  b2(t)(  J*  b2(s)ds)r  < 


for  all  r  and  t,  theorem  4.2  applies.  That  is,  if  g^It,...),  0  <  j  <  »  are  defined 


so  that 


m  m 

a(t,s)dy(s)  l  I*(k.  )  -  l  r(?. ) 

i»0  z  1  i-0  * 

/J  a(t,s)dy(s)  is  the  best  rth  order  estimate  if  and  only  if 

g^ (t,s^ ,  •  * *,Sj )  “  E{b(t)b(s1 ) •••b(s^ > >  0  <  j  <  r 

From  (3.14),  we  may  easily  calculate 


g0(t)  -  o 


g^(t,»»»)  “  j(a(t,«)  ®(t)kj_1>(  ••••) 

+  (a(t,«>  •••)  j  >  0 


(3.29) 


k  v-,-4 


v" 


•  I*  ' 


!‘iV4#>5|S3»-f  ■  ^  v-  • 


However 


j( a(t,  •)  0(t)k^1)(81,....,«j) 


”  jT  l  *(t'S(i)>Bt.n,  b<%(,)» 


l  a(t.s  )E{  n  b(s  )} 
i-1  Wi  1 


(alt,  •)  ©j(fc)k  Hs^,  •  ••,s^)  »  a(t,o)E{b(a)b(*^)...b(8_j)}do  .  (3.30) 


The  kernel  equations  (3.29)  become 


0  -  Eb(t) 


(3.31) 


alt, a)  +  /!j  a(t,o)E{b(  o)b(s)  }do  -  k{b(t)b(s)} 


(3.32) 


1  a(t,s .  )E{  n  b(s ,) }  J*  a(t,o)E{b(c)  It  b(s,)}da 


-  E{b( t)b(a^)»»*b(8j)}  ,  J  >  2 


(3.33), 


(3.31)  la  true  by  definition,  and  (3.32),  by  hypothesis.  It  remains  to  prove  that 
(3.33) j,  j  >2  all  hold.  However,  a  direct  application  of  lemma  3.5  shows  that 


E{b(o)  il  b(s  )}  »  \  mint  c,o.  )E{  It  b(s.)) 

a _ e  *  * 


for  every  j.  Using  this,  the  left  hand  aide  of  (3.33)j  becomes 
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£  {a(t,s.  )  +  f*  a(  t,  o)min(  o,s ,  )do}  S{  It  b(s  )} 

1-4  *  U  1  * 


■  J  mints  ,t)*{  It  b(s  )} 

i-1  1  l<i 

Vi 


E{b(t)b(«1)***b(«^)> 


where  the  first  aquallty  employs  the  hypothesis  on  a(t, a) ,  and  the  second  employs  lemma 
3.5  again.  Thus  (3.33).  is  true  for  all  j  >  2. 


3.4  Quadratic  Filters 

As  a  further  example  of  the  technique  of  section  3.2,  we  will  present  the  optimal 
kernel  equations  for  the  quadratic  case  (r  «  2)  and  sketch  a  theoretical  approach  to 
their  solution.  To  guarantee  validity  of  the  discussion,  assume  throughout  the  hypotheses 
of  theorem  3.2  for  r  ■  2. 

Deriving  the  optimal  kernel  equations  is  simply  a  matter  of  calculation.  Let 
f(t)  -  aQ(t)  +  jjj  •1<t,s)dy<s)  +  /Jj  Jg1  a2(t,s1,s2)dy(s2)dy(s1)  and  let  g^lt,*”)  be 
defined  from  aQ,  a 1 ,  Sj  in  the  manner  indicated  at  <3.13).  Thus 

g„(t)  -  a„<t)  +  a,(t)  k,  +  a2<t,»)  S>2  k?  (3.34) 


(3.35) 


(3.36) 


g^( t,s)  •  a^(t,s)  a0<t)k1<s)  ♦  (a^t,*)  ©1  k2>(s) 

♦  <a2<t,*)  a1  k()(s)  +  <a2(t,0  92  k3)(s) 
g2<t,sl(s2)  -  a2<t,s1,s2)  +  a(J<t)k2<81,s2)  +  (a^t,*)  9  k^ts^Sj) 
*  (a^t,*)  9^  kjXs^Sj)  +  2(a2<t,«)  ®1  ’c2><*r,2> 

+  <a2<t,*)  »2  k4)(st,s2)  . 


(More  properly,  9  in  (3.34)  -  (3.36)  should  be  5>(t).>  According  to  theorem  3.2  f(t) 
is  optimal  quadratic  iff  gQ,  g1(  and  g2,  are  respectively,  Ef(t),  Ef(t)hta)  and 
Vf(t)h{Sj  )h(s2).  Meaning  the  definition  of  9(t)  from  section  2  and  kj  » 

Eh(Sf) . • .h( Sj) ,  we  derive  for  the  optimal  kernel  equations! 


1 


t  *  *1 

Ef(t)  -  «0<t)  +  Jg  »t  <t,s)Eh(s)ds  +  /0  a2(t,»1  .s2)Eh(*1  )ht«2)d*2da1 

Ef(t)h(s)  »  at(t,s)  +  aQ(t)Eh(8)  +  af  ( t ,  o)Eh(  o)h(  s)do 


(3.37) 


*  *  °i 

+  JQ  a2(t,8,o)Eh(  o)do  +  Jg  /Q  a2(t,(J1,o2)Eh(c1)h(o2)h(8)d8 
Ef(t)h(a1)h(s2)  -  a2it,ay,a2)  +  agftJEMs,  )h(s2 )  +  a1 (t,a1 )Eh(s2 ) 


(3.38) 


+  a1(t,a2)Eh(a1)  +  J^j  a^t,  o)Eh (  o)h(81  )h(s2)do 

♦  /g  [a2(t,s1 ,  o)Eh(  o)h(s2 )  +  a2(t,a2,o)Eh(  oihla^ldo  (3.39) 

a- 

+  Jg  fQ  «2(t,a1 ,  o2)E{h(a1  )h(a2)h(  o1  )h(  o2>  }do2do1  . 

These  equations  deserve  some  elementary  remarks  before  we  set  about  solving  them. 
First,  the  optimal  kernels  are  all  interrelated  in  the  general  case.  He  cannot  solve  for 
ag  and  a^  independently  of  knowing  a2.  Likewise,  if  aQ  •  cQ,  a 1  -  c1  are  the  kernels 
of  the  beat  linear  estimate,  they  will  not,  in  general,  be  the  lower  order  kernels  of  the 
best  quadratic  estimate.  Secondly,  the  equation  (3.37)  -  (3.39)  can  be  used  for  other 
suboptimal  designs  in  the  spirit  of  theorem  3.3.  Thus,  if  aQ  and  a  ,  are  qiven,  and  we 
seek  the  best  quadratic  correction  to  «g(t)  ♦  a^ (t,s)dy( s) ,  this  will  he  found  by 

solving  (3.39)  for  a2  in  terms  of  a,  and  ag.  The  methods  developed  for  solving  the 
full  set  of  equations  will  also  apply  to  the  best  correction  problem. 

As  a  system  of  Integral  equations,  (3.37)  -  (3.39)  looks  complicated  and  contains 
unusual  features.  Nevertheless,  we  will  show  that  solving  the  system  can  be  reduced  to 
two,  familiar  tasks  —  solving  a  linear  estimation  problem  and  solving  a  Fredholm  Integral 
equation.  The  method  behind  this  reduction  is  simply  to  eliminate  aQ  and  a^  to  obtain 
an  equation  for  a2.  The  basic  steps  are:  1)  eliminate  Sg(t)  from  (3.38)  to  derive  the 
Integral  equation  (3.41)  for  a^i  2)  solve  this  for  a1  in  terms  of  a2  using  the 
solution  to  the  linear  filtering  problem,  (see  3.42)i  3)  use  (3.42)  to  eliminate  a1  from 
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(3.39)  and  derive  (3.43),  an  integral  equation  only  involving  the  unknown  a2,  and*  4) 
turn  (3.43)  into  the  Fredholm  equation  (3. 45).  The  central  equation  is  thus  (3.45).  One; 
it  is  solved  for  a2,  a,  and  aQ  are  found  by  using  (3.42)  and  (3.37)  respectively. 

Let  R  :  L^( [0,tl )  *  L^((0,tl)  be  the  operator  defined  by 


(R0)(s)  -  covth(s)  ,h(  c)l  8(c>do  . 


The  first  step  is  easy;  simply  solve  (3.37)  for  aQ(t)  and  substitute  the  result  in 
(3.38).  We  then  derive 

[I  R) a, ( t ,  •)  ( s)  -  cov(f(t),h(s)J  -  8h(o)a2(t,s,o)do 

(3.41) 

-  fQ  fQ  covfh(s),h(  ot  )h(  o2)la2(t,o1,o2)do2do1  . 

The  next  step,  solving  this  for  a.,,  thus  requires  inverting  I  +  R. 

Lemma  3.6 

i)  h(s),  s  <  t,  has  a  best  linear  estimate  h(s>  -  o^ts)  +  /*  a(s.o<  ’v(s) 

(As  a  convention,  set  a(s,o)  “  0  for  0  <  a  <  a  <  t) 
it)  (I  +  R]-1«1-Q  where  Q  is  the  Integral  operator  with  kernel 

q,sl's2)  "  «<*|t«2)  ♦  atSj.Sj)  -  alo.s^^o.SjJdo 

0  <  s1#  s2  <  t  . 

ft  2  4  — 

Proof,  we  are  assuming  E{J  h  (s)ds]  <  ®.  This  guarantees  that  h(s)  exists,  and,  as 


in  (3.25), 


a, 

a(sl(s2)  +  f  0(8^  ,c)cov(h(s2 )  ,h(  c)]d<j  ”  cov(h(s1 )  ,h(  a)) 


0  <  s2  <  s1  <  t 


11)  is  standard.  See,  for  instance,  Gessey  [3]. 

This  lemma  can  now  be  spoiled  to  solve  (3.41)  for  a.(t,s)t 


a^t/S)  «  cov(f{t),h(s)l  -  q(s,o)cov[f(t),h(o)]do 
-  %  It  .<»,«  c2)a2(t,c1,c2)da2do1 


(3.42) 


where 
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it  is  possible  to  rewrite  (3.43)  as  a  Fredholm  integral  equation  and  thereby  to  reduce  the 
task  of  calculating  a2  to  a  familiar  problem.  First  notice  that  (3.43)  may  be  rewritten 
in  the  form 


or 


«2(*1.«2>  “  F(. 


1 


>  “  <Ra2<B2' 
•  6  £  Vi 


• ) )  ( s  ^ )  -  (Ra^fSj , •) ) (Sj) 

,  s2 ,  o1 ,  a2  )*2 ( o1 ,  Oj  )dc2do1 


1(1  +  R)a2<a2> •)] (Sj)  -  F(a1 ,s2)  -  (Ra^s, , •) ) (Sj) 


"  ft  ^0  Vl'V  VWl' Wl 


(3.44) 


In  these  equations,  the  argument  t  has  been  omitted  for  simplicity.  How  apply  (X  +  R)-1 
to  both  sides  of  (3.44).  Again,  an  equation  of  the  form 

1(1  +  R]a  (st,*)l  (s2)  “  linear  terms  in  a2 

is  obtained,  but  this  time  there  are  no  partial  tensor  contractions  of  the  form  Ra^,*) 
(Sj)  on  the  right  hand  side.  With  a  final  application  (1  +  R)-1  -  I  -  Q  to  both  sides 
the  following  Fredholm  equation  for  a2  is  derived. 


jtt.s^Sj)  -  F1(t,s1,s2)  *  ft  ft  y(t,e  (,  *2'<Va2>a2,t'<V  ®2>de,2d<,1  (3.45) 


where 


|(t,s1,s2) 


F(t,Sj  ,s2)  - 
+ 


ft  lq<s2'°1  '°2>+ql81 '02>r<t'°l'82>ld02dai 

ft  ft  qC1.a1)q(s2,o2)F(t,c1,a2)dc2dc1 


y(t,B^,B2>a1,a2)  ”  r1(t,s1,s2,  <r1 , c2>  -  ft  q(s,u)Y1(u,s2,o1,o2)du 
I1(t,s1,s2,c1  ,a2)  •  -  r2(t,s1,s2,o1,c2)  -  q(s2 ,  ^  )q(s ^ ,  o2 ) 

+  io  *<Vu)r2(t,8i  (U.O^Ojjdu  . 
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Remarks .  The  viewpoint  here  is  not  recursive.  Rather  t  is  fixed  throughout  and  integral 

operators  are  defined  and  Inverted  on  L2((0,t])  or  L2{[0,tl2l,  and  at  a  later  time  t 

the  whole  operation  would  have  to  be  repeated.  This  poses  an  interesting  question  for 

further  research.  What  structure  on  the  moments  this),  Ef(t)h(s),  etc.,  would  allow  a 

recursive  solution  to  the  quadratic  kernel  equations,  in  the  sense  that  aft  *  dt,  s,,s2) 

could  be  constructed  in  a  simple  way  from  a(t,Sj,e2)?  »  related  question  is  also 

important.  When  are  the  solutions  a^  and  a2  separable  functions?  If  separability 

1  2 

occured,  then,  as  mentioned  above,  the  stochastic  Integrals  l^fa^)  and  *2^*2*  coul<J  be 
realized  as  the  outputs  of  stochastic  differential  systems.  Certainly,  if  F  and  y  of 
the  Fredholm  equation  for  a2  are  separable,  a2  will  be  separable,  but  due  to  the 
complicated  manner  in  which  the  moments  Ef(t),  Ef(t)h(s>,  etc.,  combine  to  produce  F 
and  y,  this  does  not  lead  to  easy  conditions.  This  issue  is  not  pursued  further. 
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